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ABSTRACT 

We present indications, based on a sample of 303 Abell clusters, for a relation 
between the dynamical state of clusters and the alignments of galaxy members 
with their parent cluster major axis orientation as well as with the large scale 
environment within which the clusters are embedded. The statistical results are 
complemented with a deep, wide-field case study of galaxy alignments in the 
cluster A521, which is characterised by multiple merging events (Maurogrdato et 
al 2000, Ferrari et al. 2003) and whose galaxy members show a strong alignment 
signal out to ~ 5 h^ 1 Mpc. Our results show that galaxy alignments appear to 
be stronger the more dynamically young is the cluster, especially when found in 
high-density environments. This relation complements the recently found "Clus- 
ter Substructure - Alignment Connection" (Plionis & Basilakos 2002) by which 
dynamically young clusters, found in high-density environments, show stronger 
cluster-cluster alignments. 

Subject headings: galaxies: clusters: general - large-scale structure of universe 



1. Introduction 

A well established alignment effect is that between the orientation of nearby clusters 
and between the orientation of the BCG (or cD) and that of their parent cluster (cf. Sastry 
1968; Binggeli 1982; Carter & Metcalfe 1980; Struble 1990; West 1989; West 1994; Plionis 
1994; Fuller, West & Bridges 1999; Kim et al. 2001; Chambers, Melott, Miller 2002). 
Furthermore, recent evidence has shown that substructures in clusters are aligned with the 
cluster orientation as well as along the large-scale filamentary structures within which they 
are embedded (cf. West, Jones & Forman 1995, Plionis & Basilakos 2002, Plionis 2001). 
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Analytical (cf. Bond 1987) and numerical work (cf. West, Villumsen, Dekel 1991; van 
Haarlem & Van de Weygaert 1993, Tormen 1997, Onuora & Thomas 1999; Splinter et al. 
1998; Falthebacher et al. 2002) show that cluster-cluster and substructure-cluster alignments 
occur naturally in hierarchical clustering models of structure formation, like the CDM. This 
fact has been explained as the result of an interesting property of Gaussian random fields that 
occurs for a wide range of initial conditions and which is the "cross-talk" between density 
fluctuations on different scales. This property is apparently also the cause of the observed 
filamentariness observed not only in "pancake" models but also in hierarchical models of 
structure formation; the strength of the effect, however, differs from model to model. 

Recent observational evidences point also to the so called "Cluster Alignment - Substruc- 
ture Connection" by which the dynamically young clusters show a stronger tendency to be 
aligned (Plionis & Basilakos 2002) while they are also located preferentially in high-density 
environments (Plionis & Basilakos 2002; Schiiecker et al 2001). Such alignment effects could 
be imprinted also in smaller scales; the scales of individual galaxies. For example, if galax- 
ies formed after the collapse of their parent cluster, then the anisotropic initial conditions 
could be imprinted in member galaxy orientations. In the hierarchical structure formation 
models, galaxy alignments could originate from a combination of different mechanisms, for 
example as a result of the parent cluster tidal field (cf. Salvador-Sole & Solanes 1993; Usami 
& Fujimoto 1997; but see also Barnes & Efstathiou 1987), a possibility supported also by 
the correlation found between disk galaxy spin axes and the local tidal shear field (Lee & 
Pen 2002), and/or if the galaxy-galaxy interactions occur in a preferred direction, for ex- 
ample along the primordial large-scale filamentary structure within which the protocluster 
is embedded. This anisotropic merger scenario of West (1994) has provided an interesting 
explanation of the observed strong alignment of the brightest cluster galaxy (BCG) not only 
with the cluster position angle but also with the near large-scale structures. 

In clusters, any primordial galaxy alignments effect should be severely damped by violent 
relaxation, by the exchange of angular momentum in multiple galaxy encounters that occur in 
the dense cluster environment over a Hubble time (cf. Coutts 1996) and even by secondary 
infall (Quinn & Binney 1992). Therefore, other than the alignment of the BCG with its 
parent cluster that could be expected (cf. Struble 1990; West 1994; Kim et al. 2001), it 
seems secure to say that in highly relaxed clusters, where there has been sufficient time to mix 
the phases, one should not expect to observe any significant primordial galaxy alignment, 
even if they did originally exist. Therefore, the existence or not of galaxy alignments in 
and around clusters could be an indication of their dynamical state. The possible existence 
of such alignments, intrinsic in nature, is important also for their effect on measures of 
weak lensing, where ellipticity correlations of background galaxies are expected to arise from 
the lensing of foreground large mass inhomogeneities (cf. Bartelmann & Schneider 2001). 
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Considerable effort has been devoted, recently, in attempts to quantify and disentangle these 
effects (cf. Catelan, Kamionkowski & Blandford 1996; King & Schneider 2003; Heymans & 
Heavens 2003). 

Observational efforts to detect intrinsic galaxy alignments in clusters or superclusters has 
not reached a clear consensus. For example, Dekel (1985) using the UGC and ESO-Uppsala 
galaxy catalogues concluded that no significant galaxy alignments exist in clusters. However, 
earlier studies of Adams, Storm & Storm (1980) had found alignments, although not very 
significant, in A999 and A2197. Djorgovski (1983) found galaxy alignments in the Coma 
cluster, an effect more prominent for red galaxies, while van Kampen & Rhee (1990) and 
Trevese, Cirimele & Flin (1992) analysing the 10-20 brightest galaxies in a large number of 
Abell clusters (Abell 1958; Abell, Corwin & Olowin 1989) found no galaxy alignments other 
than between the BCG and its parent cluster orientation (see Djorgovski 1987; Cabanela & 
Aldering 1998 and references therein for reviews of early results) . 

One should be very careful in how to interpret such conflicting results. It should be 
noted that systematic biases, like PSF inhomogeneities (cf. Brown et al 2002), and projection 
effects can to mask any true alignment signal. For example, the projection of foreground 
galaxies along the line of sight of a cluster as well as the projection on the plane of the 
sky of member galaxies, always work in the direction of smearing alignments. On the other 
side, telescope tracking problems can create false alignments. Such problems, if not modeled 
adequately, could potentially introduce artificial galaxy-galaxy alignments, and thus hamper 
attempts to measure weak shear induced by gravitational lensing, but cannot create a false 
alignment between galaxies and their parent cluster orientation, since the later is estimated 
from the distribution of the galaxy member positions. 

Therefore, the question whether galaxy alignments exist in clusters and if they do then 
under which circumstances do they appear, is still unanswered. In this paper we present 
a statistical analysis of galaxy alignments in 303 Abell clusters (Abell 1958; Abell, Corwin 
& Olowin 1989), based on their POSS scans and a SExtractor (Bertin & Arnouts 1996) 
identification of the galaxy distribution within a radius of 1 hr 1 Mpc from the cluster centre. 
These statistical results are complemented with a deep, wide-field case study of galaxy 
alignments within ~ 5 h^ 1 Mpc of the cluster A521, which is a relatively distant (z ~ 0.25), 
richness class 1 Abell cluster with clear evidence of multiple merging events (Maurogordato 
et al. 2000; Arnaud et al. 2000; Ferrari et al. 2003). 
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2. Galaxy Alignments in Abell Clusters 

We have analysed the POSS scans in the direction of 303 Abell clusters, all of which 
have an independent determination of their position angle (from Plionis 1994), while 168 
also have velocity dispersion measurements (Struble & Rood 1999; De Propris et al. 2001). 
Our sample, which was defined on the basis of availability of position angle measurements 
in Plionis (1994), has (z) ~ 0.07 and m w < 17. 

An area of radius spanning from ~ 17' to 40' was used so as to correspond to ~ 
1 h^ 1 Mpc at the distance of the clusters. In the cases of nearby clusters (z& 0.045) the 
corresponding radius is < 1 hr l Mpc. An interactive procedure was followed in which 
the parameters of the SExtractor identification were varied so as to maximise the number of 
detected objects, while not breaking-up the brightest galaxies into subcomponents. Although 
the SExtractor star-galaxy separation is quite efficient, it fails to correctly separate stars from 
galaxies for saturated objects and at the faint luminosity end. Bright stars were therefore 
excluded by hand. Similarly, satellite tracks and other noise signal were excluded after 
careful visual inspection of each field. 

The total number of Sextracted galaxies is ~ 5 x 10 5 out of which the majority covers 
a very small surface area and thus have unreliable classification. In the presence of no 
systematic effects, related to the size of the galaxies, we expect that the mean ellipticity 
should be independent of the galaxy image size. In figure 1 we plot the derived galaxy 
ellipticity as a function of the number of pixels per galaxy. It is evident that for 7V p i x < 30—35 
there is a systematic offset of the mean ellipticity to higher values. Therefore, in order to 
avoid biases, having to do with determining galaxy shape parameters using small number 
of pixels and the inefficiency of the star-galaxy separation at lower surface brightness, we 
choose to analyse only the higher 5% quantile of the galaxy size distribution (estimated 
independently in each field). The selection of this limit was based on a further trial and error 
approach in which the distribution of galaxy position angles was tested for homogeneity (see 
tests below) and it is such that always iVp ix > 30 — 35. 

Furthermore, we have tested, for the above size limits, the ellipticities of the unsatu- 
rated stellar images in our fields, for which deviations from zero should reflect possible PSF 
anisotropics. We find that the stellar ellipticity frequency distribution peaks at e ~ 0.05 
and has an exponential drop thereafter (figure 2). We have verified by visual inspection of 
a random subset of our fields that the major contribution to the high-e tail of the stellar 
ellipticity distribution comes from misclassified galaxies, double stars as well as background 
projections near the stellar images. 

Therefore, from the above analysis and in order to avoid ill definitions of position an- 
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gles (due to near sphericity) and possible PSF anisotropies while not throwing away too 
many galaxies, since a considerable fraction of ellipticals have relatively small ellipticities, 
we feel secure to use, in our study, galaxies that have ellipticities > 0.1. In any case, we 
expect that any residual anisotropy that still remains in our sample will be identified in 
the Fourier- analysis test, presented below, but should bare no correlation with the parent 
cluster orientations and thus should not introduce any artificial alignments. These selection 
procedures reduce the total number of analysed galaxies to 15560, with no galaxies having 
an area smaller than ~ 10" x 10" (ie., N pix > 35). 

The statistical analysis of alignments between galaxy and cluster orientations depends 
on defining the relative orientation of the galaxies with respect to the cluster, 

^ C =%-Q c \ . (1) 

where 6i and 9 C are the position angles of the galaxies and cluster, respectively. In an 
isotropic distribution we will have (0j, c ) — 45°. A significant deviation from this value would 
be an indication of an anisotropic distribution which can be quantified by (Struble & Peebles 
1985): 

5 = £^- 45 - ( 2 ) 

% 

In an isotropic distribution we will have (5) ~ 0, while the standard deviation is given 
by a = 90/Vl2iV. A negative value of 5 would indicate alignment and a positive one 
misalignment. 

2.1. Tests for systematics directional biases 

Before we analyse the possible galaxy-cluster alignments, we test whether the galaxy 
position angles, 6>j, in our final galaxy sample, have any significant directional bias. We 
quantify possible deviations from isotropy by a variety of methods: 

• estimating the Fourier transform of the galaxy position angles (cf. Struble & Peebles 
1985): 

C n = (2/iV) 1/2 cos 2nd i S n = (2/iV) 1/2 J2 sin 2n9 i ■ ( 3 ) 

If the #j's are uniformly distributed between 0° and 180°, then for large N both C n and 
S n have zero mean and unit standard deviation. This test measures the tendency of the 
position angles to cluster around some particular values. For example, the C\ and Si 
components measure deviations in the 90°, 180° and 45°, 135° directions respectively. 
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Indeed, if we find Fourier components with values > 3 this would imply systematic 
directional bias at a > 3<r level. 

• applying three statistical tests on the distribution of 0's to test whether they could 
originate from an isotropic parent distribution of position angles. These are the classical 
X 2 (parametric), the Kolmogorov-Smirnov (non-parametric) and the Rayleigh tests (cf. 
Mardia 1972). The later estimates the parameter: 

R2 = -L{Sf + C 2 ) 1/2 . (4) 

For large N we have that 2NR 2 ps \ 2 f° r 2 degrees of freedom and the probability Pr 2 
that a value of R 2 exceeds a critical value is given in Mardia (1972). 



In figure 3 we present the distribution of galaxy position angles for our total sample as 
well as for the 10 largest galaxies (excluding the BCG) in each cluster field. The distribution 
appears flat with no directional bias, as quantified also by the different statistical tests, 
presented in table 1, where we also list the values of the first Fourier moments for the 
different galaxy samples. We do find | Ci,2,3,4 1 , |<Si,2,3,4| ^ 2.3, indicating the existence of no 
significant directional bias. However, had we imposed no galaxy size limit we would have 
found significant deviations from isotropy with very large values of the first moments due 
to inefficient star-galaxy separation (we deduce this from the fact that we find the same 
anisotropy pattern when we analyse the stars in our fields). 

The above analysis is based in treating all galaxies together, irrespectively of the cluster 
field to which they belong. This does not exclude the possibility that some cluster fields 
maybe severely affected by biases, a possibility that could be hidden when treating all the 
galaxies together. To this end we calculate the Fourier moments independently in each field. 
If there is no bias then the distribution, over the 303 fields, of the individual Fourier moments 
should be a Gaussian with zero mean and unit standard deviation. In figure 4 we plot their 
distributions and the corresponding zero-mean Gaussians. We find (d) = —0.09, a Cl = 1.20 
and {Si) = 0.00, a si = 1.03. A x 2 test shows that the distributions are statistically equivalent 
to zero- mean Gaussians at a > 0.1 significance level (we find similar results for the higher 
moments as well). Therefore it seems that the small excess of Ci's at negative values is not 
very significant although it is responsible for the small negative shift of the (Ci) and for the 
20% excess of its variance. Indeed if we exclude the 7 cluster fields with \C\\ > 3 we reduce 
the overall value of C\ (see table 1) from -2.3 to -0.9 and the mean over the 296 cluster fields 
to (Ci) — —0.02. However, our alignment results, presented below, remain unchanged when 
excluding these fields. 



-7- 



2.2. Results 

In order to maximise the alignment signal to noise ratio we have stacked all cluster im- 
ages after rotating them so as to have a common orientation along their determined position 
angles. We then estimate the misalignment angle between every galaxy and the common 
cluster orientation. In order to identify where does the possible alignment signal come from, 
we apply our test to different subsamples of galaxies. One is only the Brightest Cluster 
galaxy (BCG), another is the next 10 largest galaxies (lOLGs) and finally all the galaxies 
but the BCG and lOLGs. The alignment results for different subsamples are presented in 
table 2. We find that the over all 15560 galaxies in the 303 clusters, alignment signal is 
quite weak but significant (S = —0.9, a = 0.2). Confining ourselves to the lOLGs we find a 
stronger alignment signal. Finally and as expected, we find a very strong and statistically 
significant alignment signal between the BCG's and the clusters position angles. 

2.2.1. Statistical quantification 

In figure 5 we present the distribution of misalignment angles, <fi for the BCG and 
lOLGs case. We performed a x 2 test with 5 degrees of freedom to assess whether the 
observed distribution could be drawn from an isotropic one. Furthermore we also applied a 
non-parametric Kolmogorov-Smirnov (KS) test and compared the un-binned 0-values with 
the isotropic case. The results of these tests are presented in Table 2. There are significant 
deviations from isotropy for the sample containing all the galaxies as well as for the BCG and 
the lOLGs samples. Note however, that when excluding the BCG's and lOLG's no significant 
deviation from isotropy is found (last row in table 2). For fainter galaxies the combination 
of projection effects, which tend to smear out any existing alignment, and shape parameter 
uncertainties, could hide an existing alignment effect and thus we cannot derive any strong 
conclusion from our present analysis. Unfortunately, we do not have colour information in 
order to attempt to decontaminate the cluster fields by selecting the red-sequence galaxies, 
as we will do in our case study of A521. 

Also note that the joint probability of the observing the alignments of the BCG and 
10LGS samples should increase if one considers (a) the two samples as being independent of 
each other and (b) that the deviation from isotropy, in both, is toward the same (alignment) 
direction. 
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2.2.2. Dependence on environment 

We investigate further whether there is any correlation of the alignment signal with the 
environment. It is already known that clusters in high-density environments show stronger 
cluster-cluster alignments and more substructure than in low-density ones, as would be 
expected if such clusters are still accreting matter along the large-scale filamentary structures 
in which they are embedded (cf. Plionis & Basilakos 2002; Schiiecker et al. 2001). We have 
therefore found which of the clusters in the subsample analysed do belong to superclusters. 
To this end, we have performed a percolation analysis to identify superclusters in the whole 
Abell sample with a magnitude cut m 10 < 17, corresponding to a volume limited sample 
out to z ~ 0.09. We do find a weak tendency of increasing galaxy-cluster alignment signal 
in clusters belonging to higher overdensity superclusters. For example using clusters in 
superclusters found with a percolation radius of 15 h~ x Mpc one finds (5bcg = —11.3, 
cbcg = 2.2 and 5iolg = —2.1, <tiolg = 0.7). This tendency, weak as it may be, suggests a 
correlation between the galaxy-cluster alignment signal and the cluster dynamical state, since 
in high density environments, clusters are expected to be dynamically active as discussed 
previously. 

2.2.3. Dependence on cluster dynamical state 

We want to test the alignment behaviour of the cluster high velocity dispersion tail, 
expected to correspond to highly virialised massive clusters, since the virialisation process 
should have mixed the orientation of galaxies members, and thus weak or no galaxy align- 
ments are expected. For this purpose we have analysed for alignments a subsample of 132 
clusters that have estimated a v with more than 10 galaxies per cluster. 

For the 32 clusters in our sample with a v > 900 km/sec, we obtain <5i lg — —6.1, 
(Jiolg — 1-5 with a distribution of misalignment angles being consistent with the isotropic 
case at a significance level of -P> x 2 ~ 2 x 10~ 4 and Pks ^ 10 -4 . Correspondingly for the 
remaining 100 clusters of our sample with a v < 900 km/ sec we obtain a weak alignment signal 
5 10 lg — —1-5, o iolg — 0.8. These results show a clear correlation of the alignment effect 
with velocity dispersion, which is in the opposite direction expected if the high dispersion 
clusters were indeed highly relaxed clusters. 

If we further confine our analysis in clusters belonging to dense superclusters (percolation 
radius of r perc = 15 h^ 1 Mpc) and with a > 900 km/sec (in total 19 clusters) we obtain an 
even stronger alignment signal: #iob = —7.4, Uiolg — 1-9, higher than in both cases taken 
separately. In figure 6 we present the distribution of misalignment angles for the latter case 
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(red hatched histograms) and for the clusters with a v < 900 found in any environment. 
The excess alignment signal for the high-a clusters in dense environment is evident. Their 
probability of consistency with an isotropic distribution is Pks < 0.002 for both the BCG 
and lOLGs distributions. 

At this point, the high density environment of these clusters makes us suspect that their 
high velocity dispersion does not reflect high mass virialised clusters but rather the result of 
subclustering effects, (for example, due to the large relative velocities of the different sub- 
clumps), so that these objects would be in reality dynamically young clusters. We have then 
been looking in the literature for each individual cluster if signatures of dynamical youth 
can be detected in their X-ray, optical or radio images. We have found that for 12 (A1228, 
A2061, A399/A401, A1736, A1775, A2065, A401, A426, A754, A2256 and A85) of the 19 
clusters, at least one signature does exist, for a further 3 we have not found any relevant 
information while only 4 appear to be virialised objects. 

In detail, A1228, A2061, A399/A401, A1736 and A1775 show a multi-peaked velocity 
distribution, A2065, A401, A426 and A754 show features in the X-Ray temperature maps, 
A2256 shows sub-clustering in the optical density map. Moreover, A1190, A2256 and A754 
show radio halo and relics, which are thought to be related to merging events (cf. Feretti 
& Giovanninni 1996, Feretti 2001). Finally, A85 presents signatures of merging in optical 
(presence of a foreground group in the velocity distribution), X-ray (double configuration of 
the isointensity contour map, with a primary and a secondary cluster) and in radio bands 
(presence of a radio relic). 

Since there is a triple correlation between alignments, velocity dispersion and superclus- 
tering, we present in figure 7 the alignment signal as a function of cluster velocity disper- 
sion and percolation radius (for the 132 clusters of our sample that have measured a v and 
N z > 10). The alignment signal is represented as a grey scale, with increasing brightness 
meaning more alignments, as indicated by the printed values of (5), while the contour lines 
engulf regions of some significance level (the continuous and broken lines represent the 3 and 
2.5 a levels respectively). It is evident that for the 10LG case there is a significant trend 
of increasing alignment signal both with increasing a v and supercluster percolation radius. 
It is also interesting that for the BCG's there seems to be a dependence of their alignment, 
with their parent cluster orientation, mostly to the environment, almost independent of their 
velocity dispersion. This could indicate a different (or "special") formation process for the 
BCG, with respect to the other luminous cluster galaxies (cf. West 1994). 

Our results indicate that cluster alignments are stronger in clusters located in a high 
density environment and characterised by a high-velocity dispersion. Since we already have 
independent indications that clusters in high-density environments show significantly more 
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substructure (Plionis & Basilakos 2002; Schiiecker et al. 2001), our results suggest that the 
strongest alignment effects have been detected in dynamically young clusters, ie., charac- 
terised by both a strong subclustering and a high-velocity dispersion. 

3. A case study: the dynamical active Abell 521 

The galaxy cluster A521 is a dynamically very active cluster with multiple evidences 
of merging both in optical and X-Ray data. In particular, its very high velocity dispersion 
of ~ 1300 km/s (Maurogordato et al. 2000) has been explained as the result of various 
merging events (Ferrari et al. 2003). The main component of the cluster shows a clear 
North- West /South-East direction, both in X-ray and in optical, which is the main direction 
of the undergoing merger event, but also a high density ridge in the core region perpendicular 
to the main direction, suspected to result from an older merger. The cluster is embedded 
in a dense large-scale environment (Arnaud et al. 2000). At the light of the results of the 
previous section, this high velocity dispersion and highly substructured cluster embedded in 
high density environment should be an excellent target for detecting galaxy alignments. 

Part of a larger program of a multi- wavelength follow up of merging clusters, Abell 521 
has been covered with B- and I-band imaging obtained with the CFHT12k camera over a 
field of 45x30 arcmin 2 centred on the cluster. The galaxy catalogue (positions, magnitudes 
and shape parameters) was built with SExtractor. From the analysis of the Ix(B-I) colour- 
magnitude diagram we have extracted all the galaxies lying on the red sequence region 
allowing to select the most likely early type galaxies belonging to the cluster (see Ferrari 
et al. 2003). Such a filtering of the galaxy catalog is very important since foreground or 
background galaxies can significantly smear out any true alignment effect. The alignment 
analysis presented in the present paper has been performed using the shape parameters 
measured in the B-band image leading to a more precise measurement of the position angles 
than in the I-band where the light profiles of the galaxies are more peaked. In figure 8 we show 
the stellar ellipticity frequency distribution which peaks at e ~ 0.03 and drops exponentially 
thereafter. Therefore, to avoid the gross effects of possible PSF inhomogeneities and to ensure 
a robust position angle determination (since it is ill defined for nearly spherical objects), we 
have imposed a minimum ellipticity of e = 0.1. 
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3.1. A521 Shape Parameters 

The shape parameters of the cluster (ellipticity and position angle) were estimated using 
the familiar moments of inertia method applied on the smoothed or the discrete galaxy 
distribution (cf. Basilakos, Plionis & Maddox 2000 for such an application). First, all 
galaxies within an initial small radius around the cluster cD or X-ray centre (Arnaud et al. 
2000) are used to define the initial value of the cluster shape parameters. Then, the next 
nearest galaxy is added to the initial group and the shape is recalculated and thus we obtain 
the cluster shape parameters as a function of cluster-centric distance. Within 1 h^ 1 Mpc of 
the central cD galaxy the cluster position angle and ellipticity have value of: 

# c ^140°±15° & e « 0.45 ± 0.1 , 

where the uncertainty is derived from the fluctuations around the mean by varying the 
limiting magnitude or by weighting each galaxy either by its L5 or by unity. 

3.2. Galaxy member orientations 

In figure 9 we present the galaxy position angles for different magnitude slices and for 
two limiting distances from the cluster centre, parametrised by which corresponds to 
the knee of the B-band optical galaxy luminosity function where —19.7 and therefore 

m* PS 20.8. There is an evident concentration around 120° ~ 160° in all magnitude bins, 
although at fainter magnitudes there is a slight shift toward the lower end of this range. 
This can be quantified by estimating the preferred galaxy orientation of the galaxy position 
angle distribution using the fundamental harmonics (n = 1) of their Fourier transform (cf. 
Djorgovski 1983): 

9 PI = X - arctan fa/d) . (5) 

The values of 6 pr for the different magnitude slices within 1.5 and 2.5 h~ l Mpc from the 
cluster centre, are reported in table 3 and it is evident that they are in good agreement with 
the cluster position angle which also coincides with the direction (0 ~ 138°) to the nearest 
Abell cluster neighbour (Abell 517). These results are robust to changes of the limiting 
galaxy ellipticity as can be seen in table 3. Furthermore, it is extremely interesting that the 
same galaxy preferred orientation is evident within 1.5 and also within 2.5 h^ 1 Mpc from 
the cluster central cD galaxy. The fact that the preferred galaxy orientation, 6 pr , coincides 
with the orientations of the cluster major axis, of the brightest cluster galaxies and with the 
direction to the nearest cluster neighbour is a strong indication that the anisotropic orien- 
tation distribution is not due to possible systematics and strongly supports the anisotropic 
merger structure formation scenario. 
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In table 3 we also present the outcome of the statistical tests used to assess the sig- 
nificance of our results. It is evident that the brightest and faintest slices present results 
significantly different from the Poisson expectations, while the intermediate slice shows low 
significant difference. This however could be because in this magnitude range, there appears 
to be an excess of galaxies with position angles at ~ 30°, which coincides with the direction 
of the ridge structure SI of Arnaud et al. (2000) and which is mostly evident in this mag- 
nitude range. Note that if we use all the galaxy position angles down to m ~ 24 all three 
tests show a significant difference from the Poisson expectation at a 99.999% level. 



3.3. A521 groups orientation 

Since the matter distribution outside the cluster itself could be anisotropically dis- 
tributed, not necessarily following the orientation indicated by the cluster major-axis, we 
have estimated the preferred orientation of the high-density regions as well as the galaxy 
major axis preferred orientation in independent shells of 1 Yr 1 Mpc width each. To this end 
we estimate the major axis orientation of the projected high-density excursion set (which 
correspond to groups of galaxies and different substructures) by Gaussianly smoothing the 
projected galaxy distribution on a N x N grid (with N ranging from 300 to 400) and measur- 
ing the position angles of dense groups (p/ (p) > 14) containing more than 12 cells (surface 
area Z 2700 kpc 2 ). 

We find that the group orientation distribution closely follows that of the individual 
galaxies with an almost identical preferred orientation angle, 9 pr . This is clearly seen in 
figure 10 were we present the histogram of the group and galaxies position angles (red 
and black histogram respectively), for two cluster-centric distances (one within the nominal 
radius of 2 Yr x Mpc and one between 2 and 5 h~ l Mpc; were the last group is found). Note 
that, in order to detect groups in the outer shell, we have reduced the density threshold to 
p/(p) > 11. In table 4 we present the preferred orientations of groups and of galaxies in the 
different shells, their misalignment angle and the probability that this misalignment angle 
is due to chance. The > 30° misalignment angle between the galaxy and group preferred 
orientations in the 3 rd and 4 th bin reduces to < 10° if we apply our analysis to a lower limiting 
magnitude (mbjim ~ 23.8), since we do expect a higher fraction of projection contamination 
from apparently faint galaxies outside the dense cluster region. 
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3.4. Outcome 

We therefore conclude that the dynamically very young and in the phase of merging 
Abell 521 shows significant red-sequence galaxy major-axis alignments, within the cluster but 
also in its nearby environment 5 h~ x Mpc), with (a) the cluster major-axis orientation, 
(b) the groups/substructures preferred orientation and (c) with the direction of the nearest 
Abell cluster (A517). 

4. Summary 

We have set out to investigate whether galaxies in clusters show any alignment of their 
major axis with the cluster orientation and whether such an effect is related to the cluster 
environment and dynamical state. We tackle this question by using a statistical approach, 
analysing the POSS scans of 303 Abell clusters ((z) ~ 0.07, m w < 17) and a case study of a 
highly dynamically active cluster (A521) using multicolour, wide-field, deep imaging data. 

In order not to contaminate our statistical results with the well-established BCG-cluster 
alignment we have excluded these galaxies from our analysis. We find that clusters indeed 
show significant bright galaxy major axis alignments with the cluster position angle but only 
in high-density environments (superclusters). Furthermore, the alignment signal is stronger 
in higher velocity dispersion clusters. Moreover, we have also found in our sample that the 
high velocity dispersion clusters, in high density environments, appear to be dynamically 
young, merging clusters, a fact which points strongly in the direction of galaxy alignments 
being correlated with their parent cluster dynamical youth. 

Similarly, in our detailed analysis of the dynamically very active Abell 521 we have 
found significant red-sequence galaxy alignments out to ~ 5 hr x Mpc from cluster core. The 
alignments are evident within the cluster as well in the outer environment and the preferred 
galaxy orientation is closely related to the mass distribution orientation, as well as with the 
direction of the nearest Abell cluster (A517). 

These results provide evidence that significant galaxy alignments are present in dynam- 
ically young clusters belonging in superclusters, that also extend to the large-scale environ- 
ment around clusters. Furthermore, these results lead also to the necessity of a possible 
revision of all those cluster weak lensing studies done in only one pass band, in which case 
it is impossible to separate background objects from the cluster members. 
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Fig. 1. — Galaxy ellipticities versus galaxy area (iVp ix ), used to determine shape parame- 
ters. The continuous line represents the mean ellipticity. In order to visually enhance the 
problematic region of small galaxy areas, we limit the x-axis to 100 pixels/galaxy. 
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Fig. 2. — The frequency distribution of the non saturated stellar ellipticities in our Abell 
cluster fields. 
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Fig. 3. — The distribution of galaxy position angles for the whole sample of 15560 galaxies 
(continuous line) and for the 10 largest galaxies in each cluster field 
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Fig. 4. — The distribution of the first Fourier moments of the galaxy position angles (points) 
and the Gaussian expectation in the case of no systematic directional bias. 
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Fig. 5. — The distribution of misalignment angles between the galaxy and cluster orienta- 
tions. Left panel: BCG sample. Right Panel: 10LG sample 
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Fig. 6. — The distribution of misalignment angles between the galaxy and cluster orienta- 
tions for the a < 900 km/sec clusters (black lines) and for those with a > 900 km/sec and 
found in superclusters with r perc = 15 h^ 1 Mpc (red hatched histogram). 
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Fig. 7. — Grey-scale plot of the alignment signal, (5), indicated by its value, as a function of 
membership in superclusters, of the indicated percolation radius, (y-axis) and of the cluster 
velocity dispersion (x-axis). The contours indicate the significance level at the 3 (continuous 
line) and 2.5 a (broken line) levels respectively. 
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Fig. 8. — The ellipticity frequency distribution in the A521 cluster field. 




Fig. 9. — Abell 521 galaxy position angle distribution (e > 0.1) for different magnitude limits 
within 1.5 and 2.5 h^ 1 Mpc (shaded histogram and thick red line, respectively). 
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Fig. 10. — Histogram of galaxy (black lines) and groups (red lines) position angle distribution 
in Abell 521. Left panel: The - 2 h~ l Mpc shell (m < 24.8), Right panel: The 2 - 5 /r 1 
Mpc shell (m < 23.8). 
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Table 1: Statistical tests of isotropy of the galaxy position angle distribution for the 303 
Abell cluster fields analysed. 
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Table 2: Galaxy-Cluster Alignment signal and its significance for the different samples. 
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Table 3: Abell 521 galaxy major-axes preferred orientations and alignment signal (within 1.5 
or 2.5 h^ 1 Mpc from cluster centre) in different magnitude bins. Note that the B magnitude 
completeness limit is 24 at 5a. 
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Table 4: Preferred Orientation of the Abell 521 galaxy and groups major axis distribution 
for m b < 24.8. P(< 59) is the probability that their misalignment angle is due to chance. 
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